Stability of Metallic Structure in Compressed Solid GeH4 
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We study the electronic and lattice dynamical properties of compressed solid germane in the 
pressure range up to 200 GPa with density functional theory. A stable metallic structure, Aba2, 
with a base-centered orthorhombic symmetry was found to be the lowest enthalpy phase for pressure 
from 23 to 177 GPa, suggesting an insulator to metal phase transition around 23 GPa. The Aba2 
structure is predicted to have higher superconducting transition temperature than SiH4 reported 
recently, thus presenting new possibilities for exploring high temperature superconductivity in this 
hydrogen-rich system. 



PACS numbers: 62.50.-p, 71.20.-b, 74.62. Fj, 74.70.-b 

As suggested by Ashcroft [l|, the dense hydrides of 
group IVa elements (C, Si, Ge, and Sn) are promis- 
ing candidates for realizing metallization of solid hy- 
drogen because hydrogen has already been "chemically 
precompressed" . The metallization pressures in the 
group IVa hydrides are believed to be considerably lower 
than may be necessary for solid hydrogen. Ashcroft [1] 
also suggested that these hydrogen-rich alloys are likely 
superconductors with high transition temperatures Tc. 
As put by Ginzburg [3|, "high-temperature and room- 
temperature superconductivity" and "metallic hydrogen 
and other exotic substances" are the two key "physical 
minima" at the beginning of the 21st century. Thus, 
exploring the possibility of metallic hydrogen and super- 
conductivity has long been a major driving force in high- 
pressure condensed matter science and remains an impor- 
tant challenge in modern physics. Recent experimental 
work on SiH4, using diamond-anvil cell techniques, has 
revealed an enhanced reflectivity with increasing pressure 
Q. It was found that solid SiH4 becomes opaque at 
27-30 GPa and exhibits Drude-like behavior at around 
60 GPa, signalling the onset of pressure-induced met- 
allization. After Chen's experiment, Eremets et al re- 
ported that SiH4 exhibits superconductivity at 96 and 
120 GPa[|. 

Prior to the experiment of Chen Q , several theoretical 
studies on Sill4 structure addressed the issue of whether 
the material is a favorable candidate of high tempera- 
ture superconductor. For a metallic Pman SiH4 phase, 
Feng et al [5|] obtained a superconducting transition tem- 
perature Tc of 166 K at 202 GPa by using the electron- 
phonon coupling strength for lead under ambient pres- 
sure. Pickard and Needs JJ also studied the structural 
properties of Sill4 and mentioned the possibility of super- 
conductivity in a C2/c phase. Later, Yao et al. [7| showed 
that the Pman structure is in reality not stable and that 
a new C2/c structure is dynamically stable from 65 to 
150 GPa. This C2/c Sill4 phase was predicted to exhibit 



superconductivity close to 50 K at 125 GPa 3- Recently, 
Chen et al. Q carried out first principles calculations on 
Sill4 to obtain structural, electronic, and vibrational in- 
formation, and they found out six energetically favorable 
structures in which the PI, Cmca, and C2/c structures 
have layered networks with eight-fold SiHg coordination. 
This layered feature favors metallization and were pre- 
dicted to have Tc in the range of 20 and 80 K. 

Compare to Sill4, there exist very few theoretical stud- 
ies and almost none experimental work was done on ger- 
mane. Recent progress on SiH4 make it profitable to per- 
form high pressure studies on GCH4. In early ab initio 
studies of GeH4 , Martinez et al 9] found that at around 
72 GPa the metallic SnF4-like structure becomes pre- 
ferred. Li et al [iO| followed studies of Feng et al on 
Sill4 and investigated a few possible structures. They 
found that T2 and O3 could be metallic under pressure. 
It is not known whether there exists a common structural 
feature that favors metallic Gell4, and no study on the 
possible superconductivity has been attempted. 

In this Letter we report a theoretical study of metal- 
lic phases and possible superconductivity in compressed 
solid GeH4. Considering the fact that germanium has a 
small band gap than silicon, GeH4 is expected to become 
metallic easier than Sill4. By calculating possible struc- 
tures of different crystal systems including monoclinic, 
orthogonal, tetragonal, hexagonal, and cubic structures, 
we found five energetically favorable structures in which 
the C2/c, Aba2, and Ccca structures have layered net- 
works. This layered feature favors metallization and su- 
perconductivity. The Aba2 phase is predicted to have 
considerable higher superconducting transition temper- 
ature Tc at lower pressure (<50 GPa), suggesting that 
GeH4 is another good candidate for high-temperature su- 
perconductivity. 

To study the structural and electronic behavior of 
GeH4 over a wide range of the pressure, we used gener- 
alized gradient approximation (GGA) density functional 



and projector augmented wave method as implemented 
in the Vienna ab initio simulation package (VASP) [ll|. 
An energy cutoff of 540 eV was used for the plane wave 
basis sets, and 16x16x16 and 8x8x8 Monkhorst-Pack 
k-points grids were used for Brillouin zone sampling of 
two GeH4 molecular cells and four GeH4 molecular cells, 
respectively. When searching stable structures, we per- 
formed calculations with relaxation of cell volume, cell 
shape, and ionic positions. Forces on the ions were calcu- 
lated through derivatives of the free energy with respect 
to the atomic positions, including the Harris-Foulkes like 
correction. Iterative relaxation of atomic positions was 
stopped when all forces were smaller than 0.01 eV/A. All 
possible structures were optimized using conjugate gra- 
dient scheme. The lattice dynamical and superconduct- 
ing properties were calculated by the Quantum-Espresso 
package [l^ using Vanderbilt-type ultrasoft potentials 
with a cut-off energy of 30 Ry and 300 Ry for the wave 
functions and the charge density, respectively. 24x 12 x 24 
Monkhorst-Pack fc-point grid with Gaussian smearing of 
0.05 Ry was used for the phonon calculations. Double 
fc-point grids were used for calculation of the electron- 
phonon interaction matrix element. 

The phase stability of GCH4 was systematically inves- 
tigated by means of first-principles total energy calcula- 
tions. Dozens of the selective structures were selected 
as the initial structures, which cover all the conipeting 
structures of SiH4 in previous literatures, d, [g, 0, Q It is 
worthy noting that /4i/a structure, a most stable phase 
over a wide pressure range (60-240 GPa) in SiH4, is un- 
stable in the whole pressure range we studied. Five low- 
enthalpy structures of GeH4, i.e., Aba2, Ccca, P2i/c, 
C2/c, and Fdd2, were found for pressure up to 200 GPa. 
Enthalpy differences as a function of pressure for five 
competing structures are plotted in Fig. 1. Three re- 
gions are clearly seen by our calculation: (i) 0-23 GPa, 
two monoclinic phases, Fdd2 and P2i/c dominate; (ii) 
23-177 GPa, it is Aba2, an orthorhombic phase; (iii) 177- 
200 GPa, it is monoclinic C2c. For germane, it is at 
about 23 GPa that an insulator to metal transformation 
occurs. 

Near the pressure of 2 GPa, numerous structures with 
monoclinic tetrahedral or orthorhombic symmetry, have 
very close enthalpy, constituting a "mixed" phase of 
GeH4 (see auxiliary material), which reflects the com- 
bined effect of exerted force and van der Waals force. 
Among them, two structures, one with Fdd2 symmetry 
and the other one with P2i/c, seem to be more stable. 
The two structures are very close in enthalpy as shown in 
the inset of Fig. 1. The Fdd2 phase is a face-centered or- 
thorhombic structure composed of eight four- fold coordi- 
nated tetrahedra. The eight GeH4 molecules in the unit 
cell are comparable and similarly ordered, so that the 
neighboring molecules are kept apart. With the increas- 
ing of pressure, the distance between Ge atoms decreases. 
Consequently, the shape of the whole cell changes and it 
results in the more stable P2i/c phase. The P2i/c struc- 
ture consists of four isolated covalently bonded GeH4 
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FIG. 1: (color online) The enthalpy versus pressure for five 
competitive structures of GeH4. The enthalpy of the Aha2 
phase is taken as the reference point. The inset shows en- 
thalpy difference between Fdd2 and P2i/c. 



tetrahedra with the H atom of one molecule pointing 
away from the H atoms of a neighboring molecule. Due 
to the symmetry, the arrangement of GeH4 tetrahedra 
in the Fdd2 structure is more orderly than the P2i/c 
structure. The Ge atoms in the Fdd2 phase forms quasi- 
two-dimensional sheets, while the H atoms are "fastened" 
by the Ge atoms. This makes it transform to other lay- 
ered structures more easily. The Ge-H bond lengths of 
these two structures are similar (1.53 A), leading to small 
difference in their unit volume. Subsequently, the trans- 
formation between them is very likely to occur in this 
pressure region. Because of the covalently bonded GeH4, 
it is not surprising that these two structures are insu- 
lating in this pressure range, and they were also found 
in the low-pressure Sill4. [8| The Fdd2 phase wins for 
pressure below 5 GPa, while the P2i/c phase takes over 
for pressure ranges from 5 to 23 GPa. The occurrence 
of mixed phases indicates that it is difficult to verify the 
structure of GeH4 at very low pressure (about 3 GPa) 
from experiments. 

Around 20 GPa, Fdd2 and P21/c are also energet- 
ically competitive, along with two new layered phases, 
Aha2 and Ccca. These four structures compete drasti- 
cally over the pressure range from 14 to 23 GPa for low 
enthalpy phase. The maximum enthalpy difference be- 
tween them is only 0.05 eV, indicating that solid Gell4 



undergoes a structure transition in this pressure range. 
Our calculations showed that the orthorhombic Fdd2 
structure transforms to other two orthorhombic struc- 
tures with Aba2 and Ccca symmetry at 14 GPa and 16 
GPa, respectively. Compared with Fdd2, the transition 
from monoclinic P2i/c phase to Aba2 and Ccca needs 
more pressure. P2i/c transforms to Aba2 structure at 
20 GPa and to Ccca structure at 23 GPa. It is interest- 
ing to note that Ccca structure transforms to P2i/c at 
high pressure (about 200 GPa) again. 

As shown, the most stable phase found in a wide pres- 
sure range (23 to 177 GPa) is a layered base-centered 
orthorhombic structure of space group Aba2. This is re- 
markably different from the case of SiH4, where metallic 
phases (C2/c, Cmca, PI, etc.) are higher in enthalpy than 
insulating phases for pressure lower than 200 GPa. For 
the Aba2 structure there are four formula units per unit 
cell. Four Ge atoms hold the Wyckoff 4a sites and 16 
H atoms occupy two Wyckoff 85 sites. The adjacent Ge 
layers are bridged by a pair of H atoms with Ge-H bond 
length from 1.935 A(20 GPa) to 1.673 A(175 GPa). The 
in-plane Ge-Ge bond length varies from 2.552 to 2.206 A. 
The calculated density for the Aba2 structure at 50 GPa 
is 5.397 g/cvD? . The most peculiar feature of this struc- 
ture is the exceptionally short H-H bond length about 
0.785 A from 20 GPa to 175 GPa, nearly retaining the 
same. The paired hydrogen atoms occupy the Ge layers 
with different orientation. When pressure increases, the 
orientation of paired hydrogen changes. 

A few words commenting on the Ccca phase are in or- 
der. This structure has been discussed by Yao et al in 
their study of SiH4 [3|. For GeH4, we found that Ccca 
phase is also a layered structure and is a good metal 
at the same pressure range as the Aba2 phase. The en- 
thalpy difference in Ccca and Aba2 is very small (0.11 eV 
at most) at a large pressure range, to 200 GPa. The 
Aba2 structure is about 0.053 eV per GeH4 unit lower in 
enthalpy than the Ccca structure at 20 GPa and about 
0.11 eV at 177 GPa, respectively. The Ccca structure 
contains four GeH4 units per unit cell. The hydrogen 
atoms between Ge layers connect the Ge and form 2D 
layers itself, which is parallel to Ge layers. Both Aba2 and 
Ccca structures are good candidates for metallic phases 
in this pressure region. 

Upon further compression, the metallic C2/c phase 
possesses the lowest enthalpy between 177 to 200 GPa. 
The C2/c structure transforms to the Ccca structure and 
the Aba2 structure at 154 GPa and 177 GPa, respectively. 
The C2/c phase is also a layered structure which belongs 
to monoclinic crystal system. The present C2/c structure 
is constructed by the two-dimensional Ge layers bridged 
by H atoms, and it is different from the C2/c structure 
proposed by Pickard and Needs Q when searching for the 
low-enthalpy structure of SiH4. Solid SiH4 also prefers 
to have this structure at high pressure Q. It is instruc- 
tive to observe that the obtained metallic structures with 
Aba2, Ccca and C2/c symmetry all have layered struc- 
tures. These structures are composed of distinct 2D lay- 
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FIG. 2: (color online) The energetically most favorable struc- 
tures computed for GeH4 in the pressure range 0-200 GPa. 
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FIG. 3: (color online) Electronic band structure and density 
of states for the Ada2 phase of GeH4 at the pressure of 30 
GPa. 



ers where the Ge atoms form a square-net bridged by 
H atoms. We thus obtained five energetically favorable 
structures for GeH4 at high pressure. The atomic ar- 
rangements for the four structures are shown in Fig. 2. 
Ccca structure is not given here since it has been dis- 
cussed in Sill4 by Yao et aZ.Q 

The energy band and density of states (DOS) of the 
Aba2 structure at 30 GPa pressure are shown in Fig. 3, 



as an example. Electron structure presents two promi- 
nent features. One is that a strong hybridization between 
Ge-p electrons and H-s electron appears in valence and 
conduction bands from -15 to -5 eV and from 5 to 15 eV. 
The other one is that Gc atoms expand the coordina- 
tion shell and use the 3d orbitals to bond with neighbor- 
ing H atoms. This case is also observed in SiH4 at high 
pressure. The DOS at the Fermi level (Np) for Aba2 
at 30 GPa attains a value of 1.1 (1.14 and 1.07 at 20 and 
50 GPa, respectively) states/eV/cell, exhibiting its good 
metallic property. Dominant contributions to the DOS 
near the Fermi level come from Ge-p electrons. While Ge- 
s, Ge-d, and H-s have minor contributions to the DOS 
at the Fermi level. The band structure reveals metallic 
character with large dispersion bands crossing the Fermi 
level (Ep) and a flat band in the vicinity of E^? from the 
r to the S point. The simultaneous occurrence of flat 
and steeps near the Fermi level has been suggested as 
favorable conditions for enhancing the electron pairing, 
which is essential to superconductivity. 

Further, the possibility of superconductivity for Aba2 
structure of GeH4 is discussed using the modified McMil- 
lan equation by Allen and Dynes. [13;] For comparison, 
we calculate the Tc at 30 and 50 GPa. The obtained 
electron-phonon coupling parameters A are 1.13 at 30 
GPa and 1.37 at 50 GPa, which are higher than the pre- 
dicted strong coupling value of 1.0[lJ, indicating a rather 
strong electron-phonon coupling (EPC) in GCH4. This 
strong EPC is correlated to the wide valence band and 
strong electron-electron interaction along with interband 
electron transfer. Using the Coulomb pseudopotential fi* 
of 0.1, 0.15, and 0.2, the estimated Tc are 26, 21, and 16 
K at 30 GPa and 34, 28, and 23 K at 50 GPa, respec- 
tively. Though the DOS at Fermi level decrease with 
the pressure increasing, Tc increases. This can be under- 



stood from the fact that the strengthen of EPC under 
pressure has advantage over the weaken of repulsion ef- 
fects between electrons. All obtained Tc values in GeH4 
are comparable to the Tc value recently reported in SiH4 
at 96 and 120 GPa .^4!| Most important is that no super- 
conductivity was found in Eremets et aPs experiment for 
SiH4 when the pressure is lower than 50 GPa. 

In summary, we have investigated the structural sta- 
bility of germane under pressure up to 200 GPa. The 
Fdd2 phase is confirmed to be a good candidate for the 
low-pressure insulating phase. Between 5 and 23 GPa, 
P2i/c is predicted to be the stable structure of another 
insulating phase. For a wide range of pressure, from 23 
GPa to 177 GPa, the Aba2 phase is confirmed to be the 
lowest in enthalpy. Calculations of band structure show 
that the Aba2 phase is a metal even at zero pressure. 
Thus the transition pressure from insulator to metal is 
23 GPa for germane, much lower than the theoretical 
metallization pressure of 50 GPa[i|, 60 GPa0], and 70 
GPa^S] for silane. At higher pressures (above 177 GPa), 
germane transforms into another metallic structure with 
C2/c symmetry. The estimated superconducting critical 
temperatures for the metallic Aba2 phase are compara- 
tive to those in SiH4, making GeH4 a potential supercon- 
ductor. 
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